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1 .  INTRODUCTION 


The  optical  turbulence  structure  parameter  appears  in  formulations  used  to 
estimate  the  effects  of  atmospheric  turbulence  on  imagery  and  electromagnetic 

(EM)  propagation.  For  many  optical  systems  corresponds  with  degradation  of 
performance  (Miller  and  Ricklin  1990;  Tatarski  1961;  Fried  1967).  Basic 

formulations  of  (Panofsky  1968;  Tatarski  1961;  Hill  1989;  Andreas  1988; 
Wyngaard  1973;  Wesley  1976;  Tunick  and  Rachele  1992)  include  the  gradient  of  the 
real  index  of  refraction  as  a  coefficient,  which  in  turn,  is  a  function  of  the 
gradients  of  temperature  and  moisture  (preferably  potential  temperature  and 
specific  humidity  according  to  Tatarski  (1961)).  Temperature  and  moisture 
gradients  can  be  approximated  from  sensible  and  latent  heat  flux  estimates,  and 
these  fluxes  can  be  obtained  from  radiation/energy  balance  formations. 

Numerous  radiation/energy  balance  models  appear  in  the  literature  (Angus-Leppan 
and  Brunner  1980;  Webb  1984;  Campbell  1985;  Pielke  1984;  Danard  et  al .  1984; 
Yamada  1981;  Carson  1987,  to  name  a  few)  varying  from  comparatively  simple  to 
academically  complex  and  requiring  different  amounts  and  numbers  of  inputs  and 
computer  capabilities.  In  this  report  we  present  a  model  (hereafter  known  as  the 
RT  model)  that  was  developed  for  imaging  and  EM  propagation  applications  when  a 
minimum  of  atmospheric  information  is  available.  This  report  provides  the 
concept  and  formulations  that  make  up  the  first  formalized  version  of  the  model. 
Other  reports  will  follow  as  the  model  evolves. 

The  primary  outputs  of  the  RT  model  are  estimates  of  sensible  and  latent  heat 
fluxes,  which  in  turn,  are  used  for  estimating  gradients  of  potential  temperature 
and  specific  humidity,  making  use  of  similarity  relations.  These  gradients  are 

then  used  to  estimate  the  optical  turbulence  parameter  cf,  needed  in  the  imaging 
and  propagation  equations. 

The  remainder  of  this  report  is  structured  into  the  following  sections:  Model 
Concept,  Model  Equations,  Computation  Procedure,  Sample  Results,  Summary  and 
Conclusions,  and  Literature  Cited. 

2.  MODEL  CONCEPT 

The  model  concept  is  a  fallout  from  the  operational  scenario  as  follows.  One  is 
interested  in  estimating  the  optical  turbulence  at  a  temporary  site  for 
different  times  of  one  day.  The  day  of  interest  is  known,  and  the  longitude  and 
latitude  of  the  site  are  known.  From  geological  maps  one  knows  the  soil  composi¬ 
tion;  and  from  the  meteorological  reports  one  has  an  estimate  of  soil  wetness; 
that  is,  dry,  moist,  or  wet.  However,  if  the  geological  and  meteorological 
reports  are  not  available,  one  makes  on-the-spot  judgments  by  examining  samples 
of  soil  from  the  surface  to  10  cm  below  the  surface.  In  addition  to  the  soil 
judgment,  one  also  makes  a  judgment  of  sky  conditions;  in  particular,  the  amount 
of  cloud  cover  (in  tenths  from  0  to  1),  an  estimation  of  cloud  height,  and  the 
density  of  the  clouds  (0  to  3)  .  From  the  above  information  one  proceeds  to 
estimate  the  sensible  and  latent  heat  fluxes  by  using  the  convention  and 
formulations  of  radiation  and  energy  balance  proposed  by  Carson  (1987)  as 
discussed  below. 
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2.1  Radiation  Balance  Concepts 

Carson  (1987)  states  that  the  net  radiation  flux  Rk  at  the  soil  surface  is  equal 
to  the  sum  of  the  net  short-wave  radiative  flux,  and  the  net  long-wave 

radiative  flux,  Rl(i;  that  is. 

The  short-  and  long-wave  fluxes  are  illustrated  in  figure  1. 


SHORT-MMVE  FLUXES  LONQ-MM/E  FLUXES 


(1-a)R  , 

1. 1 

e  ar* 

AIR 

SOIL 


Figure  1.  Schematic  representation  of  short-wave  and  long-wave  radiative 
fluxes  at  a  bare-soil  surface. 

In  figure  1  Rg^  is  the  downward  short-wave  radiative  flux  (including  both  the 
direct  solar  flux  and  diffuse  radiation  from  the  sky).  However,  some  of  the 
short-wave  radiation  is  reflected  at  the  earth's  surface  so  that 

■^SN  "  ~  **  ■^si  ' 

where  a  is  the  surface  short-wave  reflectivity,  commonly  called  albedo. 
However,  albedo  is  not  as  simple  as  it  appears  in  equation  (2);  it  is  a  function 
of  soil  type,  color,  vegetative  cover  and  the  elevation  angle  of  the  sun.  Values 
of  a  as  a  function  of  soil  and  vegetation  are  given  in  table  1,  mainly  from 
Danard  et  al.  (1984)  and  Campbell  (1977);  the  functional  form  of  a  relative  to 
the  sun's  position  is  given  in  section  3  (equation  (20)). 

Carson  (1987)  also  notes  that  if  is  the  downward  long-wave  radiation  and  a 
is  the  surface  absorptivity  to  long-wave  radiation,  then  the  net  incoming  flux 
from  the  atmosphere  is  (a  ■  Using  Stefan's  law,  one  can  write  that  the 

upward  flux  due  to  thermal  emission  at  the  earth's  surface  is  e  a  ,  where  is 

o  6 

the  effective  surface  temperature,  (  is  the  long-wave  emissivity  at  the  surface, 
and  CT  is  the  Stefan-Boltzman  constant.  Hence,  the  net  long-wave  radiative  flux, 
Rln,  is 

^LN  =  5  f 
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TABLE  1.  TYPICAL  VALUES  OF  ALBEDO,  a.  AND  EMISSIVITY,  e,  FOR  VARIOUS 
SURFACES* 


Surface 

a 

€ 

O'Neill  average  (D) 

0.25 

0.90 

Clay  pasture  (D) 

Dry  Clay  (G) 

Wet  Clay  (G) 

Sandy  clay 

(15%  moisture)  (J) 

0.15 

0.95 

Dry  quartz  sand  (D) 

Dry  sand  (G) 

Wet  sand  (G) 

Dry  moorland  (G) 

Wet  moorland  (G) 

Rock  (G) 

0.40 

0.80 

Ice  (G,M) 

0.65 

0.97 

Snow  (D,M) 

0.65 

New  snow  (G,DGL) 

0.85 

Old  snow  (G,DGL) 

D;  Deardorff  (1978)* ** *** 

DGL;  Danard  et  al.  (1984)** 

0.30 

G:  Geiger  (1961,  p.  29)** 

J:  Johnson  (1954,  p.  156)** 

M:  Mellor  (1977)** 


From  Campbell  (1977)*** 


Surface 

a 

Surface 

a 

Grass 

0.24 

Deciduous  woodland 

0.18 

Wheat 

0.26 

Coniferous  woodland 

0.16 

Maize 

0.22 

Swamp  forest 

0.12 

Pineapple 

0.15 

Open  water 

0.05 

Sugar  Cane 

0.15 

Dry  soil  (light  color) 

0.32 

*M.  Danard,  G.  Lyv,  and  G.  MacGillivray ,  A  Mesoscale  Bulk  Model  of  the  Atmo¬ 
spheric  Boundary  Laver.  Atmospheric  Dynamics  Corporation  report,  Victoria, 
British  Columbia,  Canada,  1984. 

**Complete  reference  information  is  contained  in  Danard,  Lyv,  and  MacGillivray. 

***G.  S.  Campbell,  An  Introduction  to  Environmental  Biophysics.  Springer  Verlag, 
New  York,  1977. 
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It  is  conunon  practice  to  combine  the  definition  of  f  with  Kxrchhoff's  law  so 
that  a  =  £.  However,  e  varies  with  soil  type,  vegetation,  and  snow  or  water 
cover  as  given  in  table  1. 

2 . 2  Surface  Energy  Flux  Balance  Concepts 

Carson  (1987)  writes  the  energy  flux  balance  at  the  soil  surface  as 

Rf,  •  H  *  L'E  *  ,  (^) 

where  is  the  net  radiative  flux,  H  is. the  turbulent  sensible-heat  flux,  L' E 
is  the  latent -heat  flux  due  to  surface  evaporation,  and  Gj  is  the  flux  of  heat 
into  the  soil.  The  convention  of  equation  (4)  is  shown  in  figure  2. 


ENERGY  FLUX  BALANCE 


R  -  H  ♦  L'E  ♦  G 

N  S 


N 

H 

AIR 

SOIL 


LE 


Figure  2.  Schematic  representation  of  the  energy  flux  balance  at  a 
bare-soll  surface. 

Equations  for  H,  L'E,  and  are  given  in  section  3.  Having  established  the  basic 
radiation/energy  balance  formulations  and  their  directional  conventions,  we  now 
discuss  the  equations  chosen  to  evaluate  them.  There  are  a  multitude  of 
relations  to  choose  from. 

3 .  EQUATIONS 

3.1  Short-Wave  Solar  Radiation 

The  formulations  we  use  to  compute  the  incoming  short-wave  solar  radiation  for 
cloudless  skies  are  patterned  after  Meyers  and  Dale  (1983)  and  Miller  and  Ricklin 
(1990).  These  are  then  augmented  with  empirical  results  by  Haurwitz  (1945)  to 
account  for  cloudy  skies. 

For  clear  skies  Meyers  and  Dale  (1983)  write 

1  =  To  T^T^cos  Z  ,  (5) 

where  I,,  is  the  extraterrestrial  flux  density  at  the  top  of  the  atmosphere  on  a 
surface  normal  to  the  incident  radiation,  Z  the  solar  zenith  angle,  and  T^  are 
the  transmission  coefficients  for  Rayleigh  scattering  (R) ,  absorption  by  perma¬ 
nent  gases  (g),  water  vapor  (w) ,  and  absorption  and  scattering  of  aerosols  (a). 
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lo  changes  thro”"hout  the  year  because  of  changes  in  the  earth-sun  distance  and 
is  adjusted  b"  using  the  equation 


4  =  1353  (W  ni-2)  {  1  +0.034 


COS 


[  2^  (n^  -  1  )  /  365  ]  1, 


(6) 


where  n'  is  the  Julian  day.  The  solar  zenith  angle  is  computed  by  using 
Z  =  cos"^  {  sin  (5)  sin  (D)  *  cos  (5)  cos  (D)  cos 

where  S  is  the  latitude,  D  the  declination  angle,  and  is  the  solar  hour 
angle.  Miller  and  Ricklin  (1990)  compute  the  solar  declination  angle  using 
formulations  by  Woolf  (1968) 

sin  D  =  sin  (23.4438)  sin  /3  , 


where  in  degrees,  is 

^  =  7  ♦  0.4087  sin  (7)  ♦  1.8724  cos  (7)  -  0.0182  sin  (27)  +  0.0083  cos  (27), 


and  7,  in  degrees,  is 


7  =  279.9348  +  d  . 


(10) 


The  angle  d  is  the  angular  fraction  of  a  year  represented  by  a  particular  date. 
The  angle  d  may  be  calculated  by 

d  =  (number  of  day  of  year  -  1)  (360/365.242)  .  (1-1) 


The  solar  hour  angle  H'  ,  in  degrees,  is  a  measure  of  the  longitudinal  distance 
from  the  sun  to  the  point  of  calculations  given  by 

//'  =  15  (T  -  M)  -  7  . 

where  T  is  the  Greenwich  mean  time  (GMT)  of  the  calculations,  M  is  the  time  in 
hours  after  midnight  of  the  passage  of  the  sun  over  the  Greenwich  meridian  or 
true  solar  noon,  and  tj  is  longitude,  counted  positive  west  of  Greenwich.  In 
terms  of  d  defined  in  equation  Ql).  M  is 

M  =  12.0  +  0.12357  sin  (d)  -  0.004289  cos  (d)  +  0.153809  sin  (2d) 

♦  0.060783  cos  (2d) 

An  empirical  equation  for  Tr  T^  in  equation  (5)  by  Kondartyev  (1969)  and  modified 
by  Atwater  and  Brown  (1974)  to  include  forward  scattering  is 

Tr  Tg  =  1.021  -  0.084  [  ro  (949  P  X  lO’^  ♦  0.051)f'^  , 

where  P  is  the  surface  pressure  (kPa)  ,  and  m  is  the  optical  air  mass  at  a 
pressure  of  101 . 3  kPa  given  by 
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ra  =  35  (1224  cos^  (Z)  +  1)'^^^,  where  Z  is  the  zenith  angle. 


(15) 


An  expression  for  computing  the  broad-band  transmission  of  water  vapor  absorption 
by  McDonald  (1960)  is 

=  1  -  0.077  (uzn)0-3  ,  (16) 


where  m  is  the  optical  air  mass  (defined  above)  and  u,  the  precipitable  water 
vapor,  is  determined  by  using  an  expression  by  Smith  (1966);  chat  is, 


u 


JTTTTy 


(17) 


where  Pj.  is  the  pressure  at  the  earth's  surface,  is  the  mixing  ratio,  g  is  the 
acceleration  due  to  gravity,  and  A  is  given  in  table  2. 

TABLE  2.  SEASONAL  AND  LATITUDINAL  MEAN  VALUES  OF  A 


Season 

Latitudinal 

Zone  Mee  Nl 

Winter 

Sorine 

Summer 

Fall 

Annual 

Averaee 

0-10 

3.37 

2.85 

2.80 

2.64 

2.91 

10-20 

2.99 

3.02 

2.70 

2.93 

2.91 

20-30 

3.60 

3.00 

2.98 

2.93 

3.12 

30-40 

3.04 

3.11 

2.92 

2.94 

3.00 

40-50 

2.70 

2.95 

2.77 

2.71 

2.78 

50-60 

2.52 

3.07 

2.67 

2.93 

2.79 

60-70 

1,76 

2.69 

2.61 

2.61 

2.41 

70-80 

1.60 

1.67 

2.24 

2.63 

2.03 

80-90 

1.11 

1.44 

1.94 

2.02 

1.62 

Northern 

Hemisphere 

average 

2.52 

2.64 

2.62 

2.70 

2.61 

A  simple  treatment  used  by  Meyers  and  Dale  (1983)  for  estimating  was  proposed 
by  Houghton  (1954) 

r,  =  X"*  ,  (1®) 

where  ra  is  the  optical  mass,  and  X  is  an  empirically  derived  constant  (on  the 
order  of  0.935). 

To  account  for  the  effect  of  clouds  we  introduce  a  transmission  coefficient  T(. 
derived  empirically  by  Haurwitz  (1945,  1946,  1948).  Haurwitz  (1945)  chose  to 
express  by 

T  =  .££  e  *""  (19) 

m 
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where  aa  and  b  are  constants  whose  values  were  determined  by  using  least  square 
techniques  applied  to  observed  clouds  of  different  percentage  sky  cover  and  of 
varying  density.  He  also  computed  the  ratio  of  insolation  with  partly  or 
completely  covered  sky  to  insolation  of  cloudless  skies  as  shown  in  table  3 . 


TABLE  3.  RATIO  OF  INSOLATION  WITH  PARTLY  OR  COMPLETELY  CLOUD -COVERED 
SKY  TO  THE  INSOLATION  WITH  CLOUDLESS  SKY  IN  PERCENT 


1  cloud  Amount 

1  Tenths 

1- 

3 

4-7 

8-9 

10 

0 

1 

2 

3 

H 

2 

3 

1 

2 

3 

0 

1 

2 

3 

H 

>.o 

103 

98 

88 

98 

93 

88 

92 

76 

68 

87 

80 

31 

18 

1.5 

100 

97 

90 

96 

89 

83 

92 

■a 

65 

88 

79 

30 

19 

2.0 

98 

96 

92 

94 

85 

78 

92 

63 

88 

79 

28 

20 

2.5 

90 

96 

94 

94 

92 

81 

74 

92 

74 

60 

89 

79 

52 

27 

20 

3.0 

85 

94 

92 

96 

90 

78 

69 

93 

74 

58 

89 

79 

50 

27 

21 

3.5 

81 

92 

90 

98 

87 

74 

65 

93 

73 

56 

90 

79 

47 

25 

21 

4.0 

77 

90 

89 

100 

84 

71 

61 

93 

72 

54 

90 

79 

46 

25 

22 

4.5 

73 

88 

87 

102 

82 

67 

58 

93 

72 

52 

91 

79 

44 

24 

23 

5.0 

70 

86 

85 

1J4 

80 

64 

54 

93 

72 

50 

91 

79 

81 

23 

23 

We  use  table  3  to  correct  our  estimates  of  insolation  for  cloud  cover. 

At  this  point  we  have  formulations  for  computing  the  incoming  short-wave  solar 
radiation  R^i ,  which  is  then  adjusted  for  albedo.  As  indicated  in  table  1, 
albedo  is  a  function  of  soil  type,  color,  moisture,  vegetation,  and  solar 
elevation  angle.  Although  the  first  four  effects  are  reflected  in  table  1,  the 
last  is  computed  by  using  an  approximation  by  Paltridge  and  Platt  (1976);  that 
is. 


q(Z)  =  Qi  +  (1  -  oi)  exp  l-k  (90°  -  Z)  ]  ,  <20) 

where  k  is  of  order  0.1,  Z  is  the  solar  zenith  angle,  and  is  the  small  zenith 
angle  value  of  albedo. 

3.2  Upward  Long-Wave  Radiation 

The  upward  long-wave  radiative  flux  is  computed  by  using  a  formulation  from 
Yamada  (1981). 


=  €g  a  r*  +  (1  -  £g) 

where  £g  is  the  surface  emissivity  (see  table  1),  and  a  is  the  Stefan- Boltzman 
constant.  To  evaluate  equation  (21)  one  requires  an  estimate  of  the  effective 
ground  temperature  Tg. 

This  is  done  as  follows;  We  assume  that  a  semiempir ical  expression  of  the 
sensible  heat  flux,  H,  given  by  Angus-Leppan  (1971)  (based  on  the  Penman  (1948) 
csmbinatinn  form)  is  adequate  as  a  first  approximation  for  clear  sky  dry  ground 
conditions.  Angus-Leppan' s  expression  is 

H  -  450  C  W  sin  ,  (22) 
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where  C  and  W  are  correction  factors  (see  figure  3  for  the  amount  of  cloud  cover 
and  ground  wetness,  respectively),  and  is  the  solar  elevation  angle 

=  90  -  Z)  .  For  clear  skies  and  dry  ground  C  -  W  -  1  in  equation  (22) . 


1.0 

iCLEAR  SKY 

1.0_ 

DRY  SURFACE 

0.8« 

1/2  SKY  COVER 

HIGH  CLOUD  (Cl) 

0.8« 

NO  vegetation 

FULL  SKY  COVER 

DRV  surface 

HIGH  CLOUD  (Cl) 

little  vegetation 

0.6_ 

1/2  SKY  COVER 

0.6« 

MOIST  surface 

MEDIUM  CLOUD  (AS.  CU) 
1/2  SKY  COVER 

LOW  CLOUD  (ST.SC.NS) 

0.4_ 

MOIST  SURFACE:  PLUS  WIND 

0.4  _ 

FULL  SKY  COVER 

0.2. 

surface  wet  or  covered  with 

FREELY-TRANSPIRING  VEGETATION 

MEDIUM  CLOUD  (AS.  CU) 

WET  OR  vegetated  SURFACE: 

0.2_ 

FULL  SKY  COVER 

LOW  CLOUD  (ST.  NS) 

0.0» 

PLUS  wind 

irrigated  field  in  midst  of  dry 

AREA;  OASIS:  PLUS  LIGHT  WIND 

0.0  _ 

FOG 

-0.2m 

irrigated  field  in  dry  AREA; 

PLUS  STRONG  WIND 

CLOUD  FACTOR  -  C 

-0.4J 

SURFACE 

WETNESS  FACTOR  -  W 

Figure  3.  Estimation  of  parameters  C  and  W. 


Having  an  estimate  for  H,  equation  (22),  we  compute  the  similarity  forms  for  the 
Monin-Obukhov  length  L  and  friction  velocity  u*  using 

-u*'  Cp  p  e,,  (23) 

■  - - 

where  Cp  is  the  specific  heat  of  air,  p  is  the  density  of  air, 

reference  level  potential  temperature,  k  is  Karraan's  constant  (0.4),  and  g  is 

acceleration  due  to  gravity. 


u*  =  k 


(24) 


where  x  =  (1  -  15  for  neutral  and  unstable  conditions,  and  x 

stable  conditions. 


1  +  5  ^  for 
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Solving  equations  (23)  and  (24)  iteratively  we  obtain  values  for  L  and  u*  for  the 
approximated  value  of  H.  Also,  from  similarity  forms  we  can  write 


6*  = 


(25) 


Hence,  we  can  compute  a  first  estimate  for  6*  ,  which  can  be  used  to  approximate 
Tg  since 


(26) 


where  y  =  (1  -  15  for  neutral  and  unstable  conditions,  and  y  =  1  +  5  ^  for 

stable  conditions. 


However,  to  evaluate  equation  (26)  we  require  an  estimate  of  (the  roughness 
length  for  temperature) ,  which  is  determined  from  our  modification  of  Verma 
(1989);  that  is. 


4  = 


f  1  +  s* 

[  1  -  B* 


(27) 


where 


f  y(z  -  d)  +  1  1  kB-i  » A  )  (28) 

[yCF-^  -  1  J® 


A  = 


+  2  tan'^.x 


(29) 


y  = 


(30a) 


X  = 


1 


(30b) 


Zg  -  roughness  length  for  momentum. 
B  -  Stanton  number, 
k  -  Karman's  constant  (0.4). 
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Values  for  kB"^  for  different  soil  surfaces  and  vegetation  are  given  in  figure  4. 
Having  an  estimate  of  Tg  for  dry  cloudless  conditions  we  compute  Rif 


D.  Reidel  Publishing  Company,  Dordrecht). 

3.3  Downward  Long-Wave  Radiation 

Gates  (1965)  offers  an  expression  for  downward  long-wave  radiation  for  clear 
skies  based  on  an  empirical  relationship  by  Swinbank  (1963) . 

Rlj  =  -17.09  +  1.195  cr  T*  ,  <31) 

where  R  is  in  milliwatts  per  centimeter^,  a  =  5.6697  X  10"’,  and  the  reference 
level  temperature  is  in  degrees  kelvin.  Paltridge  and  Platt  (1976),  in  turn, 
suggest  an  addition  to  equation  (31)  to  account  for  clouds,  giving  a  total 
expression  of 

=  -17.09  ♦  1.195  CT  r‘  +  0.3  CT  r*  (cc)  ,  (32) 

where  given  in  table  4,  is  the  emissivity  of  the  cloud  base  and  is  the 

temperature  of  the  cloud  base  in  degrees  kelvin.  Having  an  estimate  of  the  cloud 
height  one  can  approximate  ,  assuming  an  average  of  the  dry  and  moist  adiabatic 
lapse  rate. 
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TABLE  4. 

MEAN  SHORT-WAVE  PLANETARY  ALBEDO, 
EMISSIVITY,  Cgi 

«ci.  and 

LONG -WAVE  FLUX 

Cloud  Level 

Cloud  Type 

®ci 

Cci 

1 

Cirrus 

0.35 

0.3 

2 

Altocumulus -altostratus 

0.55 

0.9 

3 

Low  cloud  2 

0.60 

1.0 

4 

Low  cloud  1 

0.50 

1.0 

From  the  above  equations  we  can  compute  the  first  estimate  of  net  radiation, 
which  is  then  divided  into  sensible  heat,  latent  heat,  and  ground  heat.  For  dry, 
cloudless  conditions  we  assume  that  the  latent  heat  is  zero.  To  complete  the 
energy  balance  equation  we  also  need  a  method  for  estimating  ground  heat  flux. 

3 . 4  Ground  Heat  Flux 

In  our  initial  formulations  we  used  an  equation  by  Angus-Leppan  (1971).  Through 
trial  and  error  using  measured  flux  data  we  found  it  necessary  to  make  major 
changes  to  that  formulation,  resulting  in 


-  G-  >  (  T, 


'Sn 


)  sin  (  r  -  t„)  I  , 


(33) 


where 

G*  =  r  r,  (t„  *  2)  -  ]K,  sin  {  X  /  12  [  (t„  +  2)  -  t„  ]  }  .  (34) 

Tg„  -  ground  temperature  during  adiabatic  conditions  (approximately  1  h 
after  sunrise) 

t„  -  time  (relative  to  midnight)  of  adiabatic  conditions 
Ks 

/Cg  =  -  for  soil  type  and  wetness  (table  5,  Angus-Leppan  1971) 

2</k 

Ks  -  thermal  diffusivity 
K  -  thermal  conductivity 
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TABLE  5.  AVERAGE  THERMAL  PROPERTIES  OF  SOILS  AND  SNOW 


,  Ks 
"  2^ 


Material 

State 

Water* 

Content 

Ks 

K 

Sand 

Dry 

0.0 

0.7 

2.3 

Moist 

0.2 

4.2 

8.4 

Saturated 

0.4 

5.2 

7.4 

Clay 

Dry 

0.0 

0.6 

2.0 

Moist 

0.2 

2.8 

5.6 

Saturated 

0.4 

3.8 

5.4 

Peat 

Dry 

0.0 

■m 

0.25 

Moist 

0.4 

0.9 

Saturated 

0.8 

1.0 

Porosity'*' 

Snow 

Fresh 

0.95 

0.15 

3.0 

- 

0.8 

0.32 

1.6 

Compacted 

0.5 

1.7 

3.4 

'Volume  of  water  per  unit  volume  of  soil 
^Volume  of  pore  space  per  unit  volume 


(After  van  Wijk,  Physics  of  plant  environment.  North  Holland  1963) 

Equation  (34)  implies  that  the  ground  heat  flux  is  zero  approximately  2  h  after 
the  neutral  period  and  approximately  3  h  after  sunrise. 

We  now  introduce  formulations  by  Woolf  (1968)  that  we  used  to  approximate  the 
time  of  sunrise  and  sunset.  These  are  followed  by  our  ethod  for  computing  the 
morning  neutral  time.  The  solar  hour  angle  H'  tht.  relates  to  sunrise  and 
sunset  is 

H'  =  arcos  /  sin  A  -  sin  0  sin  D  1  (35) 

I  cos  <f>  cos  D  J 


where  ^  is  the  latitude  of  the  site,  D  is  the  solar  declination  angle  (see 
equation  (8))  and  A  is  the  sunrise/sunset  solar  elevation  angle.  At  ground  level 
Woolf  sets  A  -  -0.9  degree.  H'  in  equation  (35)  is  in  degrees,  but  can  be 
transformed  into  hours  by  dividing  by  15. 


H'C)  =  15  (T  -  W)  -  ri 


(36) 
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where  T  is  Greenwich  mean  time  (GMT)  and  M  is  the  time  of  the  sun's  passage  over 
the  meridian. 

The  solar  day  extends  from  M  -  H'  to  M  +  H' ,  and  hence  sunrise  time  is  M  -  H'  and 
sunset  time  is  M  +  H'  . 

We  also  require  the  time  after  sunrise  when  the  lower  atmosphere  is  in  a  neutral 
state,  that  is,  adiabatic,  or  when  H  —  0.  This  time  is  designated  tn  and  is 
referenced  to  midnight.  First  we  note  that  for  t  -  t^,  G  -  -G*  and  that  G* 
remains 

G-  =  [Tg  (t„  >  2)  -  TgJ  K,  sin  (  n/l2  [  (C„  +  2)  -  C„  ]  }  .  (37) 

For  very  dry  soil  L' E  ^  O  at  t„. 

Hence,  since  =  H  +  L' E  +  G  (38) 


Rn  reduces  to 

=  G  =  -G*  =  [  Tg  (t„  >  2)  -  rg„  ]  K,  sin  {  ir/l2  [  (  t„  *  2  )  -  ]  }  (39) 

where  is  evaluated  for  dry  soil. 

Since 


=  (  1  -  a  )  -R^,)  ,  (40) 

we  obtain,  by  equating  equations  (39)  and  (40), 

[^8  >  2)  -  TgJ  K,  sin  ((t„  >  2)  -  g|  =  -  (1  -  o)  R^,  *  (i?,,  -  R^,)  .  (41) 

For  neutral  conditions  we  approximate  Tg  -  Tg^  using  the  formulation 

l-sn  =  T-m  -  r.  (Z,  -  Z,)  ,  (42) 

where  T^  is  the  reference  level,  (Zj.)  is  temperature,  and  is  the  adiabatic 
lapse  rate. 

Knowing  Tg„  and  T^  one  can  compute  and  R^^  using  equations  (21)  and  (31);  that 
is , 

^Lt  =  «g  T’g  ^  (  1  -  Eg  ) 

and 

Rn  =  -17.09  +  1.195  a  T*  . 

Finally,  since  iJgj  can  be  evaluated  as  a  function  of  time,  one  can  find  the  time 
tj,  that  satisfies  equation  (41)  .  This  we  define  as  the  neutral  time  for  dry 
soil,  which  is  approximately  1  h  after  sunrise. 


(43a) 

(43b) 
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For  moist  or  saturated  soil,  the  above  procedure  Is  modified  to  allow  for  nonzero 
values  of  L'E  as  follows . 


The  latent  heat  flux  L'E  can  be  written  as 

L'E  =  ~L'  /)  u*  q*  . 


(44) 


where  u*  (equation  (24))  is  the  friction  velocity,  p  is  the  air  density,  L'  is 
the  latent  heat  of  vaporization  of  water,  and  q*  is  the  specific  humidity  scaling 
length  defined  by 


■3i  TEi  ’ 


where  k  is  Karman's  constant,  and 


for  unstable  conditions. 


=  1  ♦  5  ^  for  stable  conditions, 

Li 


^-  =  1  for  neutral  conditions. 


The  integral  of  equation  (45)  is 


9  -  Qr 


9  -  9r 


=  9* 

In 

V  - 1' 

IT 

+  1 

•  4 

- 

=  9* 

In 

-  l' 

1 

"JT 

*  1 

k  4 

I*'**  for  unstable  conditions, 


for  stable  conditions. 


and 


=  ^  In  ^  for  neutral  conditions. 


from  the  Clausius-Clapyreon  equation  we  can  write 

q  =  3 . 8  X  10-3  f  exp  | 

where  f  is  relative  humidity  and  K'  -  5.44  x  10^. 


t)}' 


(45) 


(46a) 


(46b) 


(46c) 


(47a) 


(47b) 


(47c) 


(48) 
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For  neutral  conditions  the  reference  level  and  ground  level  q  can  be  expressed 
as 


and 


(49a) 


(49b) 


Because  of  the  very  small  difference  between  „  and  Tg„,  on  the  order  of  0.02°, 
equation  (49)  can  be  differenced  to  give  an  approximation  for  equation  (47c) ; 
that  is , 


-  qrn  =  -3.8  x  10-=>  exp  j/C'  ‘  ‘  ^ 


(50) 


Using  equation  (47c) , 


-3.8  X  lo¬ 


in 


^  (^sn  -  n)  ^Xp  |  } 


(51) 


At  this  point  we  have  all  the  equations  necessary  for  computing  approximations 
for  R^,  H,  and  G  for  dry,  cloudless,  or  moist  ground  cloudy  conditions.  The  next 
section  provides  our  step-by-step  procedure  for  computing  the  sensible  and  latent 
heat  fluxes  using  the  equations  presented  above. 

4 .  COMPUTATIONAL  PROCEDURE 


In  this  section  a  flowchart  (figure  5)  of  the  computational  procedure  is  followed 
by  a  more  detailed  step-by-step  method. 

The  following  are  assumed  to  be  known  at  a  geographical  site  where  one  wishes  to 
compute  estimates  of  sensible  and  latent  heat  fluxes:  (1)  longitude  and  latitude 
of  the  site;  (2)  day  of  year  and  time  of  day;  (3)  reference  level  (2  m)  values 
of  temperature  (T^)  ,  pressure  (Pj) ,  relative  humidity  (fr),  and  horizontal 
windspeed  (V^)  ;  (4)  measure  or  judgment  of  cloud  cover  in  tenths  and  cloud 
density  on  a  scale  of  0  to  3  (3  being  most  dense);  (5)  albedo  of  the  soil, 
including  vegetative  cover;  (6)  composition  and  wetness  of  the  soil  (dry,  moist, 
saturated);  and  (7)  estimate  of  cloud  type  and  height. 


Step  1:  Compute  the  elevation  angle  =  90°  -  Z  of  the  sun  at  the  times  of 

interest  using  Woolf's  (1968)  formulations.  See  step  7. 

Step  2:  Compute  a  first  approximation  of  the  sensible  heat  flux  H  for  dry  ground 
and  cloudless  skies  using  Angus-Leppan  and  Brunner  (1980);  that  is, 


H  =  450C  W  sin  (for  C  -  W  1  in  figure  3)  . 


(52) 
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Step  3:  Compute  the  friction  velocity  u*  and  the  Monin-Obukhov  length  L. 


u*  =  V,  k 


z-d 

Zts 


(53) 


-(.-.iF. 


L  =  -u*  C  pT^/  kgH  . 


(54) 


(55) 


where  is  the  reference  level  windspeed,  k  is  Karman's  constant  (0.4),  Cp  is 
the  specific  heat  of  air,  p  is  the  air  density,  T  is  the  temperature  at 
reference  height,  g  is  acceleration  due  to  gravity,  and  is  the  roughness 
length  for  momentum,  where  from  Liu  et  al.  (1976) 


ln(2^)  =  -2.85  ■••1.19  In  Z,  , 
where  Z,  is  the  height  of  the  roughness  element. 

Step  4:  Compute  the  potential  temperature  scaling  constant  6‘  using 

r  =  u*^  /  kgL  , 

where  for  unstable  conditions,  T,  -  Q„. 


(56) 


(57) 


Step  5;  Knowing  T^,  $'  and  L,  compute  the  effective  temperature  Tg  at  height 
(RT  modification  of  Verma  (1989))  using 


T  -  T  - 


In  121^ 


-  l]  ,^r-d 

rrj  '"h  ’ 


T  17 

C  zV'" 

where  y  =  p.  -  15  .j.  for  unstable  conditions,  and 
d  -  displacement  height  -  0.7H* 


(58) 


4  = 


T5 


1  - 

r  1  >  b'Y 

1  -  B  • 

k  J  / 

(59) 


,  f7(z  -  d) 
\yU-=-dJ 


(60) 
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(61) 


+  2  tan'^  X  I  I  lo**  • 

Zg  -  roughness  length  for  momentum  over  vegetated  surface  -  0.14H*. 

H*  -  average  height  of  roughness  elements. 
kB”^  -  Stanton  number,  Verma  (1989). 

Step  6:  Compute  =  -17.09  ♦  1.195  a  T*  where  is  in  milliwatts  per 

centimeter"^  and  is  in  degrees  kelvin  (Swinbank  1963). 

Step  7:  Meyers  and  Dale  (1983)  compute  R^^  for  cloudless  skies  as 

iZst  =  I  =  lo  Tr  r,  r,  r.  cos  z  ,  (62) 

where  Iq  is  the  extraterrestrial  flux  density  at  the  top  of  the  atmosphere  on  a 
surface  normal  to  the  incident  radiation,  Z  is  the  solar  zenith  angle,  and  Ti  are 
the  transmission  coefficients  for  Rayleigh  scattering  (R) ,  absorption  by  perma¬ 
nent  gases  (g) ,  water  vapor  (w) ,  and  absorption  and  scattering  by  aerosols  (a). 


4  =  1353  (W  m2)  (i  +  0.034  cos  [  2  x  (n'  -  1)  /  365  )  )  , 


(63) 


where  n'  is  the  Julian  day,  and 

Z  =  cos'^  (sin  ^  sin  D  ♦  cos  ^  cos  D  cos  H')  , 


(64) 


where  ^  is  the  latitude,  D  is  the  declination  angle,  and  H'  is  the  solar  hour 
angle  (see  step  8)  . 


sin  D  =  sin  (23.4438)  sin  0 


(65) 


(")  =  7  +  0.4087  sin  (7)  >  1.8724  cos  (7)  -  0.0182  sin  (27) 
0.0083  cos  (27) 


(66) 


7  =  279.9348  +  d  , 


(67) 


where  d  is  the  angular  fraction  of  a  year;  that  is,  d  -  (Julian  day  -  1) 
(360/365.242)  (68) 


Tr  Tg  =  1 . 021  -  0 . 084  [  ra  949  P  x  lO'^  0 . 05l]^'^ 
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(69) 


where  P  Is  the  surface  pressure  (kPa) ,  and  m  is  the  optical  air  mass  at  a 
pressure  of  101.3  kPa 

m  =  35  (  1224  cos^  (Z)  ♦  1 
-  1  -  0.077  (um)0  3 

where  u  is  the  precipitable  water  vapor  (Smith  1966) , 

n  -  /  g  {X  *  1)  , 

where  P^  is  the  pressure  at  the  earth's  surface,  and  Ug  is  the  mixing  ratio. 
Values  of  X  are  tabulated  in  Smith  (1966). 

where  X  =  0.935  .  (73) 

Step  8:  Compute  sunrise,  local  noon,  and  sunset  time  using  Woolf's  (1968) 

formulations.  From  Woolf  the  solar  hour  angle  H'  that  relates  to  sunrise  and 
sunset  is 

H'  =  arcos  I  sin  A  -  sin  ^  sin  d1  (74) 

I  cos  (p  cos  D  J 


(70) 

(71) 

(72) 


where  0  is  the  latitude,  D  is  the  solar  declination  angle,  and  A  is  the 
sunrise/sunset  solar  elevation  angle.  At  ground  level  Woolf  sets  A  -  -0.9 
degree.  H'  is  in  degrees,  but  is  transformed  to  time  in  hours  by  dividing  by 
15.  The  solar  day  extends  from  M  -  H'  Co  M  *  H' . 

Step  9:  Compute  (1  -  a)  R,^  using  a  from  Paltridge  and  Platt  (1976). 

a  =  a'  ♦  (1  -  a')  exp  [-k(90“  -  Z)  ]  , 


where  k  - 
Step  10; 

Step  11; 


0.1,  Z  -  zenith  angle,  and  a'  is  albedo  at  high  solar  elevations 
Compute  R^f  using  Yamada  (1981) 

^LT  =  «8  7-*  ♦  (1  -  €g)  Rl*  • 

Compute  net  radiation  (Steps  1  through  15). 

^  =  (1  -  a)  Rg^  -  (Rif  -  Rii) 


(76) 


(77) 


Step  12;  Compute  using  (7^  -  (Zu_  +  d) )  ,  where  is  the 
temperature  in  degrees  kelvin  at  Z^,  7^  is  the  dry  adiabatic  lapse  rate.  andT^^ 
is  the  reference  height  temperature  at  the  time  t^  (relative  to  midnight)  of 
neutral  c^'nditions  (see  step  13). 


Step  13;  Compute  dry  value  of  G,  that  is.  Gj  using  a  formulation  by  Rachele  and 
Tunick  (see  equation  33) . 


Step  20:  Compute  for  cloudy  skies  using  Haurwitz  (1945). 

Step  21:  Compute  (1  -  a)  JJg*  . 

Step  22:  Compute  Rn  -  -17.09  +  1.195  a  T*  +  0.3  a  T*  (cc)  (86) 

where  is  in  milliwatts/centimeter^,  is  cloud  emissivity,  is  temperature 
at  cloud  base,  and  cc  is  tenths  of  cloud  cover. 

Step  23:  Compute  'Rising  the  wet  value  of  T^;  that  is, 

T*  *  (1  -  «,)  R^, 


Step  24:  Compute  R^  =  (1  -  a)  R^^  -  (^Lt  “ 

Step  25:  Compute  Tg  -  Tg„ 

Step  26:  Compute  moist  value  of  G  using  expression  in  step  13,  when  is  the 
moist  value  in  Angus-Leppan  (1971). 


Step  27  :  Compute  l'e  from  L' E  =  R^  -  G  -  H . 


(89) 


\Jr  frm/sl 

Figure  6.  The  Rachele/Tunick  function  W(Vr)  (a  modified  Angus-Leppan  and 
Brunner  plot) . 
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S.  MODEL  RESULTS 


The  results  of  this  study  are  based  on  two  sets  of  measured  data  (1/2-h  averages) 
collected  at  Davis,  California,  during  the  summer  of  1966  (Stenmark  and  Drury 
1970;  Brooks  et  al.  1968;  Morgan  et  al.  1970).  The  Davis  field  site,  a  flat,  5- 
ha  area,  at  17  m  elevation  above  sea  level,  is  located  about  2  km  west  of  the 
main  portion  of  the  University  of  California  at  the  Davis  Campus,  24  km  west  of 
Sacramento,  and  113  km  northeast  of  San  Francisco.  The  data  were  taken  during 
periods  when  the  surrounding  fields,  were  for  the  most  part,  crop  covered  and 
well  Irrigated,  giving,  in  effect,  homogenous  surface  conditions  with  respect  to 
temperature  and  moisture.  Advection  effects  were  considered  to  be  negligible. 
Profiles  of  wind,  temperature,  and  moisture  were  measured  with  transducers  at 
nine  levels  from  25  to  600  cm.  Raw  data  were  processed  to  give  1/2-h  average 
profiles . 

In  addition,  1/2-h  values  of  net  radiation  and  sensible,  latent,  and  soil  heat 
fluxes  were  available.  The  terrain  at  Davis  was  relatively  smooth  and  covered 
with  fescue  grass  (average  height  10  cm).  The  soil  was  assumed  to  be  peat  and 
was  moist.  Both  sets  of  data  were  collected  during  cloudless  sky  conditions. 

Figure  7  shows  the  1/2-h  average  values  of  reference  level  temperature,  relative 
humidity,  and  windspeed  for  the  2  days.  Figures  8,  9,  and  10  shows  the  measured 
and  model  values  of  sensible,  latent,  and  ground  heat  fluxes  for  both  days. 
Figure  11  shows  the  measured  and  model  net  radiation  for  both  days.  Figure  12 
shows  the  results  for  the  numerically  approximated  optical  turbulence  structure 
parameter  (in  the  visible)  for  both  days  using  measured  and  modeled  fluxes  for 
inputs . 
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Tme  (f^’srr) 


Figure  8. 


Measured  (solid  Line)  and  modeled  (dashed  line)  values  of  sensible 
heat  for  Davis,  CA,  data. 


Figure  9.  Measured  (solid  line)  and  modeled  (dashed  line)  values  of  latent 
heat  for  Davis,  CA,  data. 
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Figure  10.  Measured  (solid  line)  and  modeled  (dashed  line)  values  of  ground 
heat  flux  for  Davis,  CA,  data. 
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Figure  11.  Measured  (solid  line)  and  modeled  (dashed  line)  values  of  net 
radiative  flux  for  Davis,  CA,  data. 


Cn**2 


Figure  12.  The  optical  turbulence  structure  parameter,  C^,  for  visible 
wavelengths  for  Davis,  CA,  data.  Measured  fluxes  as  input 
(solid  line),  modeled  fluxes  as  input  (dashed  line). 


6.  DISCUSSION 

Figure  7  shows  that  day  7-13-66  had  both  higher  temperature  and  surface  relative 
humidity  values  than  day  6-2-66. 

A  discrepancy  appears  in  both  halves  in  figure  8  regarding  sensible  heat  flux 
estimations.  The  figure  shows  that  by  early  afternoon  the  observed  sensible  heat 
begins  to  decline  in  magnitude;  however,  an  upward  spike  appears  in  the  modeled 
results.  This  is  a  situation  where  the  observed  wind  velocity  actually  decreases 
(from  2  to  3  m/s  to  1.5  m/s);  however,  the  observed  heat  flux  does  not  respond 
by  growing  larger  in  magnitude.  Our  model,  however,  forces  a  jump  in  sensible 
heat  when  a  jump  in  windspeed  occurs. 

Since  the  estimated  latent  heat  flux  is  a  result  of  the  energy  balance  equation, 
some  discrepancy  therefore  appears  in  the  early  afternoon  in  figure  9  for  day  6- 
2-66.  For  day  7-13-66  the  net  radiation  estimates  (figure  11)  are  responsible 
for  the  discrepancies  seen  in  the  morning  hours. 

Finally  the  ground  heat  flux  estimates  seem  to  be  in  very  good  agreement  with  the 
observations,  considering  the  complexity  of  the  problem. 

7.  SUMMARY  AND  CONCLUSIONS 

The  purpose  of  this  study  was  to  determine  the  feasibility  of  structuring  a 
radiation/energy  balance  model  that  would  yield  estimates  of  sensible  and  latent 
heat  fluxes,  suitable  for  imagery  and  EM  propagation  assessments,  research,  and 
applications,  when  constrained  to  require  atmospheric  measurement  of  temperature, 
pressure,  relative  humidity,  and  wind  at  a  reference  height  (about  2  m)  only. 
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The  model  presented  satisfies  the  above  constraints;  however,  one  must  also  know 
the  day  of  the  year,  time  of  day,  longitude  and  latitude  of  the  site  of  interest, 
judgment  of  soil  type  and  moisture  (dry,  moist,  saturated),  and  cloud 
characteristics  (tenths  of  cloud  cover,  density,  cloud  type,  and  approximate 
height) . 

The  model  is  a  composite  of  formulations,  some  purely  physical  and  well  founded, 
others  semiempirical ,  and  a  few  that  are  strictly  empirical.  Our  major  contribu¬ 
tions  include  formulations  for  estimating  the  effect  of  wind  on  the  sensible  heat 
fliix;  an  expression  for  estimating  ground  heat  flux;  a  modification  of  Verma's 
expression  for  the  computation  of  the  roughness  length  of  temperature;  and,  most 
critically,  tying  together  all  the  formulations  and  establishing  a  calculation 
procedure . 

The  two  cases  presented  in  this  report  were  restricted  to  daytime  cloudless 
conditions  using  data  from  Davis,  California.  Results  are  most  encouraging  for 
determining  qualitative  estimates  of  throughout  the  day  over  fescue  grass. 

Finally,  it  is  only  fair  to  say  that  the  model  should  be  used  and  refined  for 
barren  surfaces  and  other  vegetative  cover.  A  joint  effort  with  the  U.S. 
Department  of  Agriculture  is  currently  planned. 
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Fort  Leavenworth,  KS  66027-7300 


HQDA 

DEAN-RMD/Dr.  Gomez 
Washington,  D.C.  20314 
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Conunander 
Phillips  Lab 

ATTN:  PL/LYP  (Mr.  Chisholm) 

Hanscom  AFB,  MA  01731-5000 

Director 

Atmospheric  Sciences  Division 
Geophysics  Directorate 
Phillips  Lab 

ATTN:  Dr.  Robert  A.  McClatchey 

Hanscom  AFB,  MA  01731-5000 

Raytheon  Company 

Dr.  Charles  M.  Sonnenschein 

Equipment  Division 

528  Boston  Post  Road 

Sudbury,  MA  01776 

Mail  Stop  1K9 

Director 

U.S.  Army  Materiel  Systems 
Analysis  Activity 
ATTN:  AMXSY-MP  (H.  Cohen) 

APG,  MD  21005-5071 

Commander 

U.S.  Army  Chemical  Rsch, 

Dev  &  Engr  Center 
ATTN:  SMCCR-OPA  (Ronald  Pennsyle) 

APG,  MD  21010-5423 

Commander 

U.S.  Army  Chemical  Rsch, 

Dev  &  Engr  Center 

ATTN:  SMCCR-RS  (Mr.  Joseph  Vervier) 

APG,  MD  21010-5423 

Commander 

U.S.  Array  Chemical  Rsch, 

Dev  &  Engr  Center 

ATTN:  SMCCR-MUC  (Mr.  A.  Van  De  Wal) 

APG,  MD  21010-5423 

Director 

U.S.  Army  Materiel  Systems 
Analysis  Activity 

ATTN:  AMXSY-AT  (Mr.  Fred  Campbell) 
APG,  MD  21005-5071 


Director 

U.S.  Army  Materiel  Systems 
Analysis  Activity 

ATTN;  AMXSY-CR  (Robert  N.  Marchetti) 
APG,  MD  21005-5071 

Director 

U.S.  Army  Materiel  Systems 
Analysis  Activity 

ATTN;  AMXSY-CS  (Mr.  Brad  W.  Bradley) 
APG,  MD  21005-5071 

Commander 

U.S.  Army  Laboratory  Command 
ATTN:  AMSLC-CG 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 

Commander 

Headquarters 

U.S.  Army  Laboratory  Command 
ATTN:  AMSLC-CT 
2800  Powder  Mill  Road 
Adelphi,  MD  20/83-1145 

Commander 

Harry  Diamond  Laboratories 
ATTN:  SLCIS-CO 

2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 

Director 

Harry  Diamond  T '>brvT'2t''ries 
ATTN:  SLCHD-ST-SP 

Dr.  Z.G.  Sztankay 
Adelphi,  MD  20783-1197 

AFMC/DOW 

Wright-Patterson  AFB,  OH  0334-5000 

National  Security  Agency 
ATTN:  W21  (Dr.  Longbothum) 

9800  Savage  Road 

Ft  George  G.  Meade,  MD  20755-6000 

U.  S.  Army  Space  Technology 
and  Research  Office 
ATTN;  Brenda  Brathwaite 
5321  Riggs  Road 
Gaithersburg,  MD  20882 


OIC-NAVSWC 

Technical  Library  (Code  E-232) 

Silver  Springs,  MD  20903-5000 

The  Environmental  Research 
Institute  of  MI 
ATTN;  IRTA  Library 
P.O.  Box  8618 
Ann  Arbor,  MI  48107-8618 

Commander 

U.S.  Army  Research  Office 
ATTN:  DRXRO-GS  (Dr.  W.A.  Flood) 

P.O.  Box  12211 

Research  Trianagle  Park,  NC  27709 
Dr.  Jerry  Davis 

North  Carolina  State  University 
Department  of  Marine,  Earth,  & 
Atmospheric  Sciences 
P.O.  Box  8208 
Raleigh,  NC  27650-8208 

Commander 

U.  S.  Army  CECRL 

ATTN:  CECRL- RG  (Dr.  H.  S.  Boyne) 
Hanover,  NH  03755-1290 

Commanding  Officer 
U.S.  Army  ARDEC 
ATTN;  SMCAR-IMI-I,  Bldg  59 
Dover,  NJ  07806-5000 

U.S.  Army  Communications -Electronics 
Command  Center  for  EW/RSTA 
ATTN:  AMSEL-RD-EW-SP 
Fort  Monmouth,  NJ  07703-5303 

Commander 

U.S.  Army  Communications  -  Electronics 
Command 

ATTN:  AMSEL-EW-D  (File  Copy) 

Fort  Monmouth,  NJ  07703-5303 

Headquarters 

U.S.  Army  Communications -Electronics 
Command 

ATTN:  AMSEL-EW-MD 

Fort  Monmouth,  NJ  07703-5303 


Commander 

U.S.  Army  Satellite  Comm  Agency 

ATTN:  DRCPM-SC-3 

Fort  Monmouth,  NJ  07703-5303 

Director 

EW/RSTA  Center 

ATTN:  AMSEL-EU-DR 

Fort  Monmouth,  NJ  07703-5303 

USACECOM 

Center  for  EW/RSTA 
ATTN;  AMSEL-RD-EW-SP 
Fort  Monmouth,  NJ  07703-5303 

6585th  TG  (AFSC) 

ATTN:  RX  (CPT  Stein) 

Holloman  AFB,  NM  88330 

Department  of  the  Air  Force 
OL/A  2nd  Weather  Squadron  (MAC) 
Holloman  AFB,  NM  88330-5000 

PL/WE 

Kirtland  AFB,  NM  87118-6008 
Director 

U.S.  Army  TRADOC  Analysis  Command 
ATTN:  ATRC-WSS-R 

White  Sands  Missile  Range,  NM  88002 
Rome  Laboratory 

ATTN;  Technical  Library  RL/DOVL 
Griffiss  AFB,  NY  13441-5700 

Department  of  the  Air  Force 
7th  Squadron 
APO,  NY  09403 

AWS 

USAREUR/AEAWX 
APO,  NY  09403-5000 

AF  Wright  Aeronautical  Laboratories 
Avionics  Laboratory 
ATTN:  AFWAL/AARI  (Dr.  V.  Chimelis) 
Wright- Patterson  AFB,  OH  45433 

Commander 

U.S.  Army  Field  Artillery  School 
ATTN:  ATSF-F-FD  (Mr.  Gullion) 

Fort  Sill,  OK  73503-5600 
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Conunandant 

U.S.  Army  Field  Artillery  School 
ATTN;  ATSF-TSM-TA 
Mr.  Charles  Taylor 
Fort  Sill,  OK  73503-5600 

Commander 

Naval  Air  Development  Center 
ATTN:  A1  Salik  (Code  5012) 
Warminister,  PA  13974 

Commander 

U.S.  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-DA-M 

Mr.  Paul  Carlson 
Dugway,  UT  84022 

Commander 

U.S.  Army  Dugway  Proving  Ground 
ATTN;  STEDP-MT-DA-L 
Dugway,  UT  84022 

Commander 

U.S.  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-M  (Mr.  Bowers) 
Dugway,  UT  84022-5000 

Defense  Technical  Information  Center 
ATTN;  DTIC-FDAC 
Cameron  Station 
Alexandria,  VA  22314 

Commanding  Officer 
U.S.  Army  Foreign  Science  & 
Technology  Center 
ATTN :  CM 

220  7th  Street,  NE 
Charlottesville,  VA  22901-5396 

Naval  Surface  Weapons  Center 
Code  G63 

Dahlgren,  VA  22448-5000 


Commander  and  Director 

U.S.  Army  Corps  of  Engineers 

Engineer  Topographies  Laboratory 

ATTN:  ETL-GS-LB 

Fort  Belvoir,  VA  22060 

TAC/DOWP 

Langley  AFB,  VA  23665-5524 

U.S.  Army  Topo  Engineering  Center 

ATTN;  CETEC-ZC 

Fort  Belvoir,  VA  22060-5546 

Commander 
Logistics  Center 
ATTN;  ATCL-CE 
Fort  Lee,  VA  23801-6000 

Commander 

USATRADOC 

ATTN;  ATCD-FA 

Fort  Monroe,  VA  23651-5170 

Science  and  Technology 
101  Research  Drive 
Hampton,  VA  23666-1340 

Commander 

U.S.  Army  Nuclear  &  Cml  Agency 
ATTN;  MONA-ZB  Bldg  2073 
Springfield,  VA  22150-3198 


Commander 
U.S.  Array  OEC 
ATTN:  CSTE-EFS 

Park  Center  IV 
4501  Ford  Ave 

Alexandria,  VA  22302-1458 
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